Near-term tropical cyclone risk and coupled Earth system model biases.

Most current climate models predict that the equatorial Pacific will evolve under greenhouse gas-induced warming to a more El Niño-like state over the next several decades, with a reduced zonal sea surface temperature gradient and weakened atmospheric Walker circulation. Yet, observations over the last 50 y show the opposite trend, toward a more La Niña-like state. Recent research provides evidence that the discrepancy cannot be dismissed as due to internal variability but rather that the models are incorrectly simulating the equatorial Pacific response to greenhouse gas warming. This implies that projections of regional tropical cyclone activity may be incorrect as well, perhaps even in the direction of change, in ways that can be understood by analogy to historical El Niño and La Niña events: North Pacific tropical cyclone projections will be too active, North Atlantic ones not active enough, for example. Other perils, including severe convective storms and droughts, will also be projected erroneously. While it can be argued that these errors are transient, such that the models' responses to greenhouse gases may be correct in equilibrium, the transient response is relevant for climate adaptation in the next several decades. Given the urgency of understanding regional patterns of climate risk in the near term, it would be desirable to develop projections that represent a broader range of possible future tropical Pacific warming scenarios-including some in which recent historical trends continue-even if such projections cannot currently be produced using existing coupled earth system models.

A mutual information criterion with applications to canonical correlation analysis and graphical models.

This paper derives a criterion for deciding conditional independence that is consistent with small-sample corrections of Akaike's information criterion but is easier to apply to such problems as selecting variables in canonical correlation analysis and selecting graphical models. The criterion reduces to mutual information when the assumed distribution equals the true distribution; hence, it is called mutual information criterion (MIC). Although small-sample Kullback-Leibler criteria for these selection problems have been proposed previously, some of which are not widely known, MIC is strikingly more direct to derive and apply.

Standard assessments of climate forecast skill can be misleading.

Assessments of climate forecast skill depend on choices made by the assessor. In this perspective, we use forecasts of the El Niño-Southern-Oscillation to outline the impact of bias-correction on skill. Many assessments of skill from hindcasts (past forecasts) are probably overestimates of attainable forecast skill because the hindcasts are informed by observations over the period assessed that would not be available to real forecasts. Differences between hindcast and forecast skill result from changes in model biases from the period used to form forecast anomalies to the period over which the forecast is made. The relative skill rankings of models can change between hindcast and forecast systems because different models have different changes in bias across periods.

Deterministic skill of ENSO predictions from the North American Multimodel Ensemble.

Hindcasts and real-time predictions of the east-central tropical Pacific sea surface temperature (SST) from the North American Multimodel Ensemble (NMME) system are verified for 1982-2015. Skill is examined using two deterministic verification measures: mean squared error skill score (MSESS) and anomaly correlation. Verification of eight individual models shows somewhat differing skills among them, with some models consistently producing more successful predictions than others. The skill levels of MME predictions are approximately the same as the two best performing individual models, and sometimes exceed both of them. A decomposition of the MSESS indicates the presence of calibration errors in some of the models. In particular, the amplitudes of some model predictions are too high when predictability is limited by the northern spring ENSO predictability barrier and/or when the interannual variability of the SST is near its seasonal minimum. The skill of the NMME system is compared to that of the MME from the IRI/CPC ENSO prediction plume, both for a comparable hindcast period and also for a set of real-time predictions spanning 2002-2011. Comparisons are made both between the MME predictions of each model group, and between the average of the skills of the respective individual models in each group. Acknowledging a hindcast versus real-time inconcsistency in the 2002-2012 skill comparison, the skill of the NMME is slightly higher than that of the prediction plume models in all cases. This result reflects well on the NMME system, with its large total ensemble size and opportunity for possible complementary contributions to skill.

Assessing probabilistic predictions of ENSO phase and intensity from the North American Multimodel Ensemble.

Here we examine the skill of three, five, and seven-category monthly ENSO probability forecasts (1982-2015) from single and multi-model ensemble integrations of the North American Multimodel Ensemble (NMME) project. Three-category forecasts are typical and provide probabilities for the ENSO phase (El Niño, La Niña or neutral). Additional forecast categories indicate the likelihood of ENSO conditions being weak, moderate or strong. The level of skill observed for differing numbers of forecast categories can help to determine the appropriate degree of forecast precision. However, the dependence of the skill score itself on the number of forecast categories must be taken into account. For reliable forecasts with same quality, the ranked probability skill score (RPSS) is fairly insensitive to the number of categories, while the logarithmic skill score (LSS) is an information measure and increases as categories are added. The ignorance skill score decreases to zero as forecast categories are added, regardless of skill level. For all models, forecast formats and skill scores, the northern spring predictability barrier explains much of the dependence of skill on target month and forecast lead. RPSS values for monthly ENSO forecasts show little dependence on the number of categories. However, the LSS of multimodel ensemble forecasts with five and seven categories show statistically significant advantages over the three-category forecasts for the targets and leads that are least affected by the spring predictability barrier. These findings indicate that current prediction systems are capable of providing more detailed probabilistic forecasts of ENSO phase and amplitude than are typically provided.

Pacific Hurricane Landfalls on Mexico and SST.

A statistical model of Northeast Pacific tropical cyclones (TCs) is developed and used to estimate hurricane landfall rates along the coast of Mexico. Mean annual landfall rates for 1971-2014 are compared to mean rates for the extremely high Northeast Pacific sea-surface temperature (SST) of 2015. Over the full coast, the mean rate and 5%-95% uncertainty for Saffir-Simpson category one and higher TCs (category-1+ TCs) is 1.241.051.33yr-1 for 1971-2014 and 1.690.892.08yr-1 for 2015, a difference that is not significant. However, the increase for the most intense landfalls, category-5 TCs, is significant: 0.0090.0060.011yr-1 for 1971-2014 and 0.0310.0160.036yr-1 for 2015. The SST impact on the category-5 TC landfall rate is largest on the northern Mexican coast. The increased landfall rates for category-5 TCs is consistent with independent analysis showing that SST has its greatest impact on the formation rates of the most intense Northeast Pacific tropical cyclones. Landfall rates on Hawaii ( 0.0330.0190.045yr-1 for category-1+ TCs and 0.0100.0050.016yr-1 for category-3+ TCs for 1971-2014) show increases in the best estimates for 2015 conditions, but the changes are insignificant according to our tests.

Long-lead prediction of the 2015 fire and haze episode in Indonesia.

We conducted a case study of NCEP CFSv2 seasonal model forecast performance over Indonesia in predicting the dry conditions in 2015 that led to severe fire, in comparison to the non-El Niño dry season conditions of 2016. Forecasts of the Drought Code (DC) component of Indonesia's Fire Danger Rating System were examined across the entire Equatorial Asia region and for the primary burning regions within it. Our results show that early warning lead times of high observed DC in September and October 2015 varied considerably for different regions. High DC over Southern Kalimantan and Southern New Guinea were predicted with 180-day lead times, whereas Southern Sumatra had lead times of up to only 60 days, which we attribute to the absence in the forecasts of an eastward decrease in Indian Ocean SSTs. This case study provides the starting point for longer-term evaluation of seasonal fire danger rating forecasts over Indonesia.

More tornadoes in the most extreme U.S. tornado outbreaks.

Tornadoes and severe thunderstorms kill people and damage property every year. Estimated U.S. insured losses due to severe thunderstorms in the first half of 2016 were $8.5 billion (US). The largest U.S. effects of tornadoes result from tornado outbreaks, which are sequences of tornadoes that occur in close succession. Here, using extreme value analysis, we find that the frequency of U.S. outbreaks with many tornadoes is increasing and that it is increasing faster for more extreme outbreaks. We model this behavior by extreme value distributions with parameters that are linear functions of time or of some indicators of multidecadal climatic variability. Extreme meteorological environments associated with severe thunderstorms show consistent upward trends, but the trends do not resemble those currently expected to result from global warming.

Human influence on tropical cyclone intensity.

Recent assessments agree that tropical cyclone intensity should increase as the climate warms. Less agreement exists on the detection of recent historical trends in tropical cyclone intensity. We interpret future and recent historical trends by using the theory of potential intensity, which predicts the maximum intensity achievable by a tropical cyclone in a given local environment. Although greenhouse gas-driven warming increases potential intensity, climate model simulations suggest that aerosol cooling has largely canceled that effect over the historical record. Large natural variability complicates analysis of trends, as do poleward shifts in the latitude of maximum intensity. In the absence of strong reductions in greenhouse gas emissions, future greenhouse gas forcing of potential intensity will increasingly dominate over aerosol forcing, leading to substantially larger increases in tropical cyclone intensities.

Tornado outbreak variability follows Taylor's power law of fluctuation scaling and increases dramatically with severity.

Tornadoes cause loss of life and damage to property each year in the United States and around the world. The largest impacts come from 'outbreaks' consisting of multiple tornadoes closely spaced in time. Here we find an upward trend in the annual mean number of tornadoes per US tornado outbreak for the period 1954-2014. Moreover, the variance of this quantity is increasing more than four times as fast as the mean. The mean and variance of the number of tornadoes per outbreak vary according to Taylor's power law of fluctuation scaling (TL), with parameters that are consistent with multiplicative growth. Tornado-related atmospheric proxies show similar power-law scaling and multiplicative growth. Path-length-integrated tornado outbreak intensity also follows TL, but with parameters consistent with sampling variability. The observed TL power-law scaling of outbreak severity means that extreme outbreaks are more frequent than would be expected if mean and variance were independent or linearly related.

Rapid intensification and the bimodal distribution of tropical cyclone intensity.

The severity of a tropical cyclone (TC) is often summarized by its lifetime maximum intensity (LMI), and the climatological LMI distribution is a fundamental feature of the climate system. The distinctive bimodality of the LMI distribution means that major storms (LMI >96 kt) are not very rare compared with less intense storms. Rapid intensification (RI) is the dramatic strengthening of a TC in a short time, and is notoriously difficult to forecast or simulate. Here we show that the bimodality of the LMI distribution reflects two types of storms: those that undergo RI during their lifetime (RI storms) and those that do not (non-RI storms). The vast majority (79%) of major storms are RI storms. Few non-RI storms (6%) become major storms. While the importance of RI has been recognized in weather forecasting, our results demonstrate that RI also plays a crucial role in the TC climatology.